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1. Motivations 



^ c"| Nonequilibrium quantum fields play a central role in various areas of 

Q^, physics, such as cosmology, ultra-relativistic heavy-ion collisions, or con- 

O !* densed matter physics. A timely example in cosmology concerns the issue 

of reheating 2 of the early universe: At the end of inflation the universe is 
left in a frozen albeit very high-energy state, very different from the subse- 
quent hot universe observed from the cosmic microwave background. This 
amounts to an enormous entropy growth due to particle production. It is 
important progress that questions such as the far-from-equilibrium particle 
production at the end of the inflationary universe, including the subsequent 
process of thermalization, can be addressed in quantum field theory from 
first principles. 3 The late-time approach to quantum thermal equilibrium 
represents a complex collective phenomenon, which can occur on time scales 
dramatically larger than typical relaxation times. 4,5 ' 6 For phenomenological 
applications it is therefore often crucial that different quantities effectively 
thermalize on different time scales and a complete thermalization of all 
quantities may not be necessary. For instance, an approximately time- 
independent equation of state p = p(e), characterized by an almost fixed 
relation between pressure p and energy density e, can form very early - even 
though the system is still far from equilibrium. Such prethermalized quan- 
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titics approximately take on their final thermal values already at a time 
when the occupation numbers of individual momentum modes still show 
strong deviations from the late-time Bose-Einstein or Fermi-Dirac distri- 
bution. Prcthcrmalization is a universal far-from-equilibrium phenomenon 
which occurs on time scales dramatically shorter than the thermal equi- 
libration time. 7 In order to establish such a behavior it is crucial to be 
able to compare between the time scales of prethermalization and thermal 
equilibration. Approaches based on small deviations from equilibrium, or 
on a sufficient homogeneity in time underlying kinetic descriptions, arc not 
applicable in this case to describe the required "link" between the early and 
the late-time behavior. 

Similar issues arise in the study of high-energy nuclear collisions, where 
highly excited matter is produced far from equilibrium. The experiments 
seem to indicate early thermalization whereas the present theoretical un- 
derstanding of QCD suggests a much longer thermal equilibration time, 8 
thereby questioning for instance the picture of a weakly interacting partonic 
gas. 9 This also points to the need to go beyond kinetic approaches, which 
do not provide a valid description on very short time scales. Efforts to go 
beyond kinetic descriptions include studies of the non-equilibrium dynam- 
ics of a chiral quark-meson model, 5,7 or photon production 10 in heavy ion 
collisions. 

The recent progress in the theoretical description of nonequilibrium 
quantum fields is based on efficient functional integral techniques, so-called 
n-particlc irreducible effective actions, for which powerful nonperturbative 
approximation schemes are available. 11 The techniques provide an efficient 
approximation scheme also for situations close to, or in equilibrium where 
resummations are needed. This includes e.g. the calculation of transport 
coefficients 12 and thermodynamic quantities 13 in hot field theories or the 
study of critical phenomena near second-order phase transitions. 14 These 
are promising systematic approaches to cure the convergence problems of 
perturbative calculations of thermodynamic quantities at high temperature 
such as pressure or entropy. 15 

Though we focus on particle physics and cosmology applications, we 
emphasize that these techniques can be equally applied to other nonequi- 
librium phenomena in complex many body systems. A prominent example 
is the study of Bosc-Einstcin condensates. The interpretation of various re- 
cent experiments, for instance concerning Bosc condensation of molecules, 16 
indicate important contributions beyond standard mean-field type approx- 
imations for the real-time dynamics of correlations. 17 These can be system- 
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atically taken into account using nPI effective action techniques. 18 

Apart from these applications, there has been much progress in under- 
standing formal aspects of approximation schemes based on nPI effective 
actions, such as renormalization and the validity of Ward identities for 
gauge theories. In the following we review these developments. 

2. Non-equilibrium quantum field theory: methods 

2.1. Standard approximations fail out of equilibrium 

Nonequilibrium dynamics requires the specification of an initial state at 
some time t = 0, characterized by an initial density matrix pu- All the 
information about the time-evolution of the system is encoded in the cor- 
relation functions of field operators ip, which may be obtained from the 
generalized generating functional: 



with n-point classical sources R n = R n (x\, . . . , x n ). Here, R n ■ (p n = 
ixi...x n ^nfaii ■■■■> x n) <p(xi) ■ ■ ■ <p(x n ), with J x = J c dx° J d 3 x. The closed 
time contour C, 19 runs from the initial time t = 0, where the system has 
been prepared, to an arbitrarily large time, and backward to the initial time. 
The functional (1) can be given the following path-integral representation: 



where S[ip] is the classical action and Dip represents the integral over 
classical paths along the contour C, with fixed boundaries on both ends at 
t = corresponding to the field configurations <po( x ) an d ( / ? o( x )- One n_ 
nally performs the appropriate average over possible initial conditions. All 
correlation functions can be given similar expressions by functional differen- 
tiation. A direct evaluation of the corresponding path integrals in real time 
is currently not possible with standard Monte-Carlo techniques. Equiva- 
lcntly, the solution of the real-time dynamics of the functional Schrodingcr 
equation is prohibitively difficult. A frequently employed strategy is to 
concentrate on classical statistical field theory instead, which can be sim- 
ulated. This can give important insights when the number of field quanta 
per mode is sufficiently large such that quantum fluctuations are suppressed 
compared to statistical fluctuations. However, classical Rayleigh-Jeans di- 
vergences and the lack of genuine quantum effects - such as the approach 




(1) 
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to quantum thermal equilibrium characterized by Bose-Einstein or Fcrmi- 
Dirac statistics - limit their use. A coherent understanding of the time 
evolution in quantum field theory is required - a program which has made 
substantial progress in recent years with the development of powerful the- 
oretical techniques. 

There exist various powerful approximation methods for vacuum field 
theory as well as for thermal equilibrium at nonzero temperature. How- 
ever, most standard field theoretical techniques are not applicable for out- 
of-equilibrium situations, where one encounters additional complications. 
The first new aspect concerns secularity: The perturbative time evolution 
suffers from the presence of spurious, so-called secular terms, which grow 
with time and invalidate the expansion even in the presence of a weak cou- 
pling. Here it is important to note that the very same problem appears 
as well for nonperturbative approximation schemes such as standard 1/N 
expansions, where N denotes the number of field components. 51 Typically, 
secularity is a not a very difficult problem and for a given approximation 
there can be various ways to resolve it. There is a requirement, however, 
which poses very strong restrictions on the possible approximations: Uni- 
versality, i.e. the insensitivity of the late-time behavior to the details of the 
initial conditions. If thermal equilibrium is approached then the late-time 
result is universal in the sense that it becomes uniquely determined by the 
energy density and further conserved charges. To implement the neces- 
sary nonlinear dynamics is demanding. Both requirements of a non-secular 
and universal behavior can indeed be fulfilled using efficient functional in- 
tegral techniques: so-called n-particle irreducible effective actions, where 
the simplest example is given by the two-particle-irreducible (2PI) effective 
action, 20 which is described in the following. 

2.1.1. Secularity 

Here, we discuss the well-known secularity problem in the context of a 
relativistic real scalar field theory with potential m 2 (p 2 /2 + A</? 4 /4!. For the 
argument it is irrelevant that we assume spatial homogeneity and isotropy 
as well as Z 2 -symmetry, so that Tr{/j£> <f(x)} = 0, with Gaussian initial 
density matrix p£>. 



a Note that restrictions to mean-field type approximations such as leading-order large- 
are insufficient. They suffer from the presence of an infinite number of spurious conserved 
quantities and are known to fail to describe thcrmalization. Secularity enters the required 
next-to-leading order corrections and beyond. 
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The time-ordered propagator along the contour C is defined as: 
G(x, y) = Tr{ PD T c (<p(x) <p(y))} = 6+ G> (x, y) + 6^ G<(x, y) , (3) 
where ef = e c (±(x° - y )). One obtains to 0(X 2 ): 

G 2 -\oo V {x,y) = G (x - y) + / Go(x-u)Y. 2 (u-v)Go(v~y), (4) 

J uv 

where the bare inverse propagator is given by Gq (x, y) = i(D+m 2 )5c(x—y) 
and where S 2 denotes the self-energy to two-loop order: Y^ 2 (x — y) = 
—i\Go(x,x)6c(x — y)/2 — X 2 Go(x — y)/6. The first (one-loop) contribu- 
tion corresponds to a simple mass shift and plays no role for the present 
argument. We ignore it as well as all local insertions in the following dis- 
cussion. We therefore write: 

V2(x-y) = ~G 3 {x-y). (5) 

For illustrational purposes, we introduce an effective particle number 
n p (t), defined as: 

w p (2np(t) + 1) = [d t d t ,G<(t,t';p)+LJ 2 p G<(t,t';p)] , (6) 



where u> p = \J p 2 + m 2 and where G(t,t';p) denotes the spatial Fourier 
transform of the propagator. Writing explicitly the time integrals in (4), 
one easily obtains, for the two-loop contribution to the particle number: 

n 2 p -^(t) - n» = (1 + nW)K<(t) - n<°>K>(t) , (7) 

where, introducing the Fourier transform in frequency space Tlf' > (ijj,p): 

*<■>(*) = -i- / % '-^-^ • (8) 

LUp J Z7T (LO — LO v Y 

In the large-time limit, one gets: 

n<>>(t)~a<'> xt, (9) 

with the on-shell rates: ct< > = — S^ :> (cj p ,p)/2u;p. We recover the usual 
secular behavior of perturbation theory which occurs whenever the rates 
in (9) do not vanish on-shcll. For instance, in the present two-loop approx- 
imation, the production rate a p receives contributions from in-medium 



two-body scatterings q + k — ► r + p: 
A 2 



- f dqdkdr (2tt) 4 S^(Q + K - R - P) x n° q n° k (1 + n°) . (10) 
'p J 

where we introduced four-momenta as e.g. = (u> q , q) and dq = j^Pj^H: ■ 
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A similar analysis reveals that higher-order contributions to Eq. (4) 
are secular as well. The standard solution to this problem corresponds to 
summing the infinite series of self-energy insertions: 

G2 = Go + Go • S2 • Go + Go • £2 ' Go • £2 ' Go + • • • 

= G + Go-£ 2 -G2, (11) 

where we have used an obvious schematic notation. This corresponds to 
the usual resummation in terms of one-particle- irreducible (1PI) proper ver- 
tices. In the two-loop example discussed here, this leads to a non-secular 
time-evolution. However, the problem reappears for higher order contribu- 
tions to the self-energy. For instance, the following 0(A 4 ) contribution 

^G 2 (x -y)x J G (x - u) G 3 (u - v) G (v - y) (12) 

contains the same time-integral as analyzed above and thus diverges at 
large times. Therefore, the systematic resolution of the secularity prob- 
lem requires infinite resummations beyond the standard 1PI resummation 
scheme. 

To further elaborate the discussion, we consider the sum of contributions 
(5) and (12), which can be written as: S 4 = -(A 2 /6)G§ {G + G • £ 2 • G }. 
Our previous analysis suggests to resum the infinite geometric series of self- 
energy (£2) insertions in the term in brackets: — (A 2 /6) G 2 , {Go + Go • £2 • 
Go + . . .} = — (A 2 /6) Gq G2. Similarly, one can resum self-energy insertions 
on the other two lines of the original two-loop diagram, to arrive at the 
following simple - and suggestive - expression: 

A 2 

£[2] = — §" G 2- (I 3 ) 

To improve on this, we note that the simple expression (13), suggests the 
following iterative procedure to include higher orders without encounter- 
ing new secular terms: From the propagator at the n-th iteration G[2„], 
compute the resummcd two-loop self-energy Spn] = — ^ 2 G[2n]/6. The next 
iteration G[2( n +i)] is then computed from G^ n+1 ^ = Gq 1 — £[ 2ll ]. Start- 
ing from G[ ] = Go, this converges, after infinitely many iterations, to the 
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following equation for G = G^: 13 

G- 1 = Go 1 - E[G] (14) 

with 

£[G] = -yG 3 . (15) 

We arrive at a non-linear equation for the two-point function G. The pro- 
cedure described above guarantees that possible secular terms are automat- 
ically resummed at each step of the iteration. Equation (15) is remarkably 
simple and suggests to reorganize the original perturbative expansion in 
terms of the full propagator. Before we discuss below that this is efficiently 
described in terms of the 2PI effective action, we describe in the following 
that this procedure also automatically respects late-time universality. 



2.1.2. Universality 

The price to be paid for the above non-secular scheme is that at some order 
one necessarily has to solve nonlinear equations without further approxima- 
tions. However, as we will see next, it is precisely the nonlinearity which is 
required to be able to obtain universality in the sense mentioned above. We 
consider first Eq. (11), which is more conveniently rewritten for initial-value 
problems as: Gq 1 • G2 = 1 + £2 • Gi. Using this equation, evaluated along 
the contour C, one obtains for the time-derivative of the effective particle 
number n p (t): 

to p n p (t) -Rc^ dt' {z<(t,t';p){id t -LJ p )G%{t',t;p) 

-X>(t,t';p)(id t -uj p )G<(t',t;p)}. (16) 

To illustrate the content of this equation, we employ the following free-field 
like ansatz: (id t — uj p )G£ (f ,t;p) — > —n p (t) e -4 "'^* ~ l \ and similarly for G^ 
with n p (t) replaced by 1 + n p (t). Using (8)) one finds: 

h p (t) = (l + n p (t))K<(t)-n p (t)K>(t). (17) 



b This procedure is a direct generalization of the successive summations of so-called daisy 
and super-daisy diagrams, which correspond to the standard Hartrec approximation. 
With present notations, the latter corresponds to iterating the one-loop self-energy: 
E[ n j (a:, a;') = — iXGy^ (x, x) &c( x — x ')/2- I n the discussion presented here, we have 
discarded these (trivial) one-loop tadpole insertions for simplicity. 
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Using the late-time behavior (9), one finally gets: 

n p (t)~n° p e-^ t + ^-(l-c-^ t ) (18) 

where 7 P = a> — a p . The late-time limit, indicated on the RHS, is com- 
pletely determined by the initial particle number n p , see (10). This depen- 
dence on initial conditions is rooted in the fact that Eq. (11) is a linear 
equation for the propagator. Late-time universality - and the associated 
effective loss of information - requires nonlinear equations. 

To illustrate the importance of the resummation procedure discussed 
in Sec. 2.1.1 and, in particular, of the non-linear character of the resulting 
equations, we consider the time evolution of the particle number (6) for 
the approximation (15). We emphasize that these particle numbers are 
only employed for illustrational purposes. The quantum field theoretical 
description for the propagator contains off-shell as well as memory effects 
and is not limited to late-times. 4 Moreover, the functional formulation of 
this approach, based on the 2PI effective action described below provides 
a systematic framework for higher-order calculations. Using Eq. (14), it 
is easy to check that, in the present approximation, the time derivative of 
n p (t) has a similar form as in (16) with the replacements G^ :> — > G <> 
and ^ ~* S < ' > . Employing a similar free-field like ansatz as above, one 
obtains the following Boltzmann-like equation at sufficiently large times: 

\2 f 

"p(*) = T~ \ dqdkdr (2tt) 4 S (4) {Q + K-R-P) 

x {n,(t)n fc (t)(l + n r {t)){l + n p (t)) 
~n p {t)n r {t){\ + n k (t))(l + n,(t))} . (19) 

The late-time equilibrium solution - corresponding to h p (t) — - is given 
by the Bose-Einstein distribution: 

^= ^J )t T_ V (20) 

and is only characterized by the temperature T and chemical potential \i. 
The latter are uniquely determined by the total energy and particle number 
densities, which are conserved in the present approximation. As expected, 



c While the former is an exactly conserved quantity, the latter is only conserved in the 
approximation (19), which describes on-shcll two-body elastic scatterings. In particular, 
it is not conserved by the original equation (14), which includes off-shell, particle-number- 
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the late-time result is insensitive to the details of the initial condition. 
Note that, in contrast, resummation schemes based on local insertions, 
such as e.g. the so-called 2PPI loop expansion, fail to describe late-time 
thermalization at finite loop order. 21 This demonstrates that it is crucial 
to go beyond local (mass) resummations. 



T[tj>, G] = S[<j>] + - Tr In G" 1 + - Tr G^ 1 • G + T 2 [<j>, G] , (21) 



2.2. The 2PI effective action 

The previous analysis illustrates that a systematic way of respecting both 
non-secularity as well as universality is to reorganize the pcrturbative ex- 
pansion in terms of the full propagator. This can be done in an efficient 
way by performing a double Legendre transform of the functional (f ) with 
respect to the linear and bilinear sources Ri(x) and R,2(x,y). One then 
obtains the so-called two-particle-irreducible (2P1) effective action, a func- 
tional of the one and two-point functions <fi(x) and G{x,y), which can be 
conveniently parametrized as: 20 

-THnG- 1 + - 
2 2 

where T2 can be written as an infinite series of closed 2P1 diagrams with 
lines corresponding to G and vertices given by the shifted action S[(f> + <p]. 
Physical solutions corresponds to the stationarity conditions: 

Q^! = „ , M =0 , (22) 

oifi dG 

Using (21), the second of these equations is equivalent to the equation of 
motion (14) for G, with self-energy given by: 

^,G]=2*^M. (23) 

Thus, one automatically obtains non-linear, self-consistent equations, as re- 
quired to describe late-time thermalization. This generalizes the resumma- 
tion procedure described previously. d Systematic approximations include a 
loop or a 1/iV-expansion of the 2PI effective action beyond leading order. 22 
Finally, the effective action T[4>], which encodes all the n-point functions 
of the theory, is obtained as: 

T[ct>} = T[cf>,G[4>}}, (24) 



changing processes. When the latter are taken into account, the equilibrium distribution 
is given by (20) with a vanishing chemical potential. 4 

d For instance, the resummed two-loop self-energy (15) is obtained from the three-loop 
contribution to the functional T2: i(A 2 /48) J xy G 4 (x,y). 
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where G[<f>] is the solution of the second of Eqs. (22). Equation (24) is 
an exact relation between the 1PI and 2PI functionals which merely states 
that both are equivalent representations of the theory in the absence of 
higher-than-linear sources R n >2- However, systematic (e.g. loop or l/N) 
expansions of the 1PI and 2PI functionals do not coincide order by order in 
general. Equation (24) therefore defines a non-trivial approximation scheme 
for the computation of n-point functions. We call it the 2PI resummation 
scheme. It is non-perturbative in nature and, as we have argued previously, 
it includes all the ingredients needed to overcome both the secularity and 
universality problems of non-equilibrium quantum field theory. 



2.3. Symmetries 

Linear symmetries of the classical action are directly promoted at the level 
of the 2PI effective action. For example, for a (3(iV)-symmetric scalar field 
theory T[cj>, G] is invariant under the transformation: 

<P(x)^K<P(x) ■ G(x,y) ^TZG(x,y)n f , (25) 

where TZ denotes a O(N) rotation. This ensures that the resummed effec- 
tive action (24) has the same symmetry properties as the classical action: 
T[(j)} = r[TZ(j)]. Notice that this holds for any 2PI approximation which 
respect the symmetry. In turn, this directly implies that Ward identities 
are automatically satisfied by the n-point functions calculated from the re- 
summed effective action (24). For instance, in the case of spontaneously 
broken symmetry, the two-point function, evaluated at the solution <f> = (j> 
of the first of Eqs. (22): 



ab{ ,V) 5<t> a {x)5<t> b {y) 



(26) 



where a, b are 0(N) indices, exactly satisfies Goldstone's theorem . 22 ' 24 

We emphasize that similar arguments can be applied to arbitrary linear 
and/or affine symmetry transformations including local ones. For instance, 
for QED in linear gauges, one can show that the 2PI resummed effective 
action satisfies the usual Ward identities. In particular, one immediately 



e In contrast, the solution G a t,[4>] of the second of Eqs. (22) violates Goldstone's theorem 
in general. However, it can be shown that for systematic expansion schemes, such as a 
2PI loop or 1 /TV-expansion at a given order, the functions (26) and iG~^[4>] only differ 
by higher-order terms and, therefore, these violations are of higher order as well. 23 
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concludes that the photon polarization tensor calculated from it is trans- 
verse in momentum space at any order of a 2PI loop-expansion. 23 The 
generalization to nonabelian gauge theories is technically more involved 
and needs to be further investigated. 



2.4. Renormalization 

A systematic renormalization procedure of the resummed propagator at 
vanishing field, G[(f> = 0], has been proposed in Refs. 24 ' 25 for the 2PI 
loop-expansion in scalar field theories/ This is based on applying a BPHZ 
subtraction procedure to diagrams with resummed propagators. In par- 
ticular, it has been shown that the subtleties arising from self-consistent 
resummations reduce in this case to a proper renormalization of the follow- 
ing Bcthe-Salpeter like equation (here, G = G[<p = 0]): 



V(x, y, z, t) = A(x, y,z,t)+ % - I A(x, y, u, v)G{u 1 u)G{v, v)V(u, v, z, t) , 

^ J UVUV 

(27) 

where 

Hx,y,z,t) = A f^=J; G \ . (28) 

5G(x,y)SG(z,t) g=g[^>=o] 

This analysis has been recently generalized in 13,23 to include the field- 
dependence of the 2PI resummed effective action (24). This allows one to 
determine all the counterterms of the theory at a given level of approxima- 
tion from suitable renormalization conditions and to obtain finite results 
for arbitrary n-point functions. As an explicit example, we consider the 
one-component scalar field theory discussed in the previous sections at or- 
der \ 2 R in the 2PI coupling-expansion, where Xr denotes the renormalized 
coupling. We employ renormalization conditions for the two-point function, 
r' 2 ), and four-point function, , given by: 

r' 2 'M),^ , (29) 

6<j>{x)6<j>(y) 0=0 

r(4)(z,y ' z,t) E ^o' (30) 

Without loss of generality we use renormalization conditions for <f> = 



f For a recent application to the 2PI 1/iV-expansion, sec 
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which in Fourier space read: 

ZT^(p 2 )\ p=0 = -m 2 R , (31) 

^r(% 2 )| F o = -l, (32) 

,P2,P3)\pi=p 2 =P3=0 — — Afl , (33) 



with the wave function renormalization Z. Defining the renormalizcd fields 
as: (f>R = Z~ x l 2 (\> and Gr[4>r\ = Z~ X G\<\?\, and introducing the usual coun- 
tcrtcrms: Zm 2 = m 2 R + 5m 2 , Z 2 \ = Xr + SX and 5Z = Z — 1, the classical 
contribution to the effective action (24) reads: 

f ( 1 + 5Z „ , „„ , m?> + <5m 2 , 9 A B + 5X lA \ 
S[<f>] = J [—j—d^aWfa - ^ _ — ^4 j _ (34) 

Similarly, one can write for the one-loop part: Tr In G -1 [<j>] = Tr In G^ 1 [4>r] 
up to an irrelevant constant, and: 

1 -TtG^[4>)G[4>] = ~ J [(1 + tfZip* + m 2 , + <5m 2 ] G R (x, y; <f> R )\ x=y 
_ X R + (5Ai y fy^Q^^.^^ ^ 

Here, (5.Zi, (5m 2 and (5Ai denote the same counterterms as above, however 
approximated to the given order. Finally, the contribution from the 2PI 
functional T 2 in Eq. (21) reads, at order X 2 R : 

V2 



J XV 



L ^ ■' xy 



One first has to calculate the solution Gr[4ir] of the stationarity con- 
dition for the 2PI effective action. For this, one has to impose the same 
renormalization conditions as for the propagator (31)-(32) in Fourier space: 

iGr 1 ^ 2 ] 4> R )\ p =o,<t> R =o = -m 2 R , (37) 

-^ % G R \ P 2 - ^)| p=0 ,^ = o = -1 , (38) 

together with the condition: 

Vr(P1,P2,P3)Ipi=P2=P3=0 = -Aij, (39) 

for the renormalizcd "four-point" field, obtained from Eq. (27). It is easy 
to check that conditions (37)- (39) fix all the counterterms of the effective 



February 2, 2008 1:14 WSPC/Trim Size: 9in x 6in for Proceedings 



paper 



13 

action at vanishing field, namely 8Z\, 8m\ and SX 2 . This corresponds to 
the case considered in 24 - 25 . 

The remaining counterterms are to be determined from the renormal- 
ization conditions (31)-(33). A similar analysis to the one performed in 
25 shows that the coupling divergences appearing in the two-point function 
(29) can be absorbed in a proper renormalization of the following four-point 
field: 23 

V R {x,y^ ^}'f^\ , (40) 
0(pR{z)(j)R{t) 4> R =o 

which satisfies a Bcthe-Salpeter like equation similar to (27). The corre- 
sponding condition reads, in momentum space: 

Vr(P1,P2,P3)Ipi= P 2=P3=0 = ~^R- ( 41 ) 

The latter, together with Eqs. (31)-(32) fix the counterterms associated 
with the quadratic field-dependence in Eqs. (34)-(36), namely 6Z, 5m 2 and 
<5Ai. g Finally, the remaining counterterm SX, associated with the quartic 
field-dependence, is determined from Eq. (33). 

It is important to emphasize that Eqs. (37) and (38), as well as (39) 
and (41) are not independent conditions. Instead, they are fixed by the 
renormalization conditions (31)-(33). This express the fact that in the exact 

(2) 

theory, the following relations between rcnormalizcd quantities hold: T R = 
iG] i 1 ((j)R = 0) and = Vr = Vr. Thus, the renormalization procedure 
described here provides an efficient fixing of all the above counterterms, 
including those associated with the explicit field-dependence of the effective 
action. The latter play a crucial role e.g. in the broken phase and are 
essential for the determination of the effective potential. 



2.5. Equivalence hierarchy for nPI effective actions 

The 2PI resummation scheme provides a self-consistent dressing of the field 
<f> and the propagator G. This can be generalized to higher n-point func- 
tions as well: The nPI effective action T[(j), G, V3, . . . , V n ] provides a self- 
consistent description for the dressed proper vertices V3, . . . , V n , which are 
obtained the stationarity conditions: ST/SV^ = 0, . . . , 5T/SV n = 0. As 



s In the present approximation, one has 8Z1 = SZ and Sm^ = 8m , which follow from 
the identity F^ 2 '(x,y) = iG~~ 1 (x,y; <j> = 0), as well as 8X2 = 8X1, which follows from 
V(x, y; z, t) = V(x, y, z, t). 13 Notice that, although these identities are true in the exact 
theory, they may not be satisfied in general at the level of approximations. 
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before, non-secular approximation respecting late-time universality can be 
obtained from systematic expansion schemes of higher effective actions. In 
general this leads to non-linear integro-differential equations for the n-point 
functions, which are now independent self-adjusting variables. As described 
previously, this is a crucial ingredient for late-time universality. 

The use of nPI effective actions with n > 2 can be important to describe 
non-Gaussian initial density matrix for non-equilibrium systems. They 
are also relevant in the context of high-temperature gauge theories for a 
quantitative description of transport coefficients beyond the "leading-log" 
approximation, 27 or in the context of critical phenomena near second-order 
phase transitions. For instance, the quantitative description of the criti- 
cal exponents of scalar (f> 4 theory goes beyond any finite order 2PI loop- 
expansion 11 . 14 It requires taking into account vertex corrections that start 
with the 4PI effective action to four-loop order. 

For practical purposes, a crucial observation is that there exists an 
equivalence hierarchy between nPI effective actions. For instance, in the 
case of a loop expansion one has: 27 

r (Hoop) ^ ^ p(lloop) ^ G ] = r (lloop) ^ g,V s ]=... , 

p(21oop) ^ ^ p(21oop) G] = r (21oop) [0 } G> y 3 ] = r (21oop) Q ,V 3 ,V 4 } = . . . 
pOloop) ^ ^ p(31oop) Qj jL r (31oop) [0 } G> y 3 ] = pOloop) [^G,V 3 ,V 4 } = . . . 

where r(" loop ) denotes the approximation of the respective effective action 
to n-loop order in the absence of external sources. E.g. for a two-loop ap- 
proximation all nPI descriptions with n > 2 are equivalent and the 2PI 
effective action captures already the complete answer for the self-consistent 
description up to this order. In contrast, a self-consistently complete re- 
sult to three-loop order requires at least the 3PI effective action in general, 
etc. 1 Typically the 2PI, 3PI or maybe the 4PI effective action captures 
already the complete answer for the self-consistent description to the de- 
sired/computationally feasible order of approximation. 

We thank the organizers for this stimulating meeting. We also thank 
Sz. Borsanyi, U. Reinosa and C. Wetterich for fruitful collaborations. 



h Critical phenomena can be described using the 2PI l/./V-cxpansion beyond LO. 14 
'There exist further simplifications which can decrease the hierarchy even further such 
that lower nPI effective actions can be sufficient in practice. For instance, for a vanishing 
field 4>, one has, in the absence of sources: r( 31oo P) [<j> = 0, G] = r( 31oo P) [</> = 0, G, V 3 = 

o, vy = . . .. 
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